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Abstract. There exist two central measures of turbulent 
mixing in turbulent stratified fluids that are both caused by 
molecular diffusion: 1) the dissipation rate D{APE) of 
available potential energy APE; 2) the turbulent rate of 
change Wr,turbuient of background gravitational potential 
energy GPEr- So far, these two quantities have often been 
regarded as the same energy conversion, namely the irre- 
versible conversion of APE into GPEr, owing to the well 
known exact equality 13 (APE') = Wr,turbuient foraBoussi- 
nesq fluid with a linear equation of state. Recently, how- 
ever, Tailleux (2009) pointed out that the above equality no 
longer holds for a thermally-stratified compressible, with the 
ratio ^ = Wr,turbuientlD{APE) being generally lower than 
unity and sometimes even negative for water or seawater, 
and argued that D{APE) and Wr, turbulent actually repre- 
sent two distinct types of energy conversion, respectively the 
dissipation of APE into one particular subcomponent of in- 
ternal energy called the 'dead' internal energy IEq, and the 
conversion between GPEr and a different subcomponent 
of internal energy called 'exergy' lE^xergy In this paper, 
the behaviour of the ratio ^ is examined for different strat- 
ifications having all the same buoyancy frequency N ver- 
tical profile, but different vertical profiles of the parameter 
T = aP/{pCp), where a is the thermal expansion coeffi- 
cient, P the hydrostatic pressure, p the density, and Gp the 
specific heat capacity at constant pressure, the equation of 
state being that for seawater for different particular constant 
values of salinity. It is found that ^ and Wr,turbuient depend 
critically on the sign and magnitude of dT/dz, in contrast 
with D{APE), which appears largely unaffected by the lat- 
ter. These results have important consequences for how the 
mixing efficiency should be defined and measured in prac- 
tice, which are discussed. 
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1 Introduction 

As is well known, turbulent diffusive mixing is a physical 
process that it is crucially important to parameterise well in 
numerical ocean models in order to achieve realistic simula- 
tions of the water mass proper ties and of the so-called merid- 



ional overturning circulation dGregg L 119870 . which are two 



essential components of the large-scale ocean circulation that 
may interact with Earth climate. For this reason, much effort 
has been devoted over the past decades toward understanding 
the physics of turbulent mixing in stratified fluids, one impor- 
tant goal being the design of physically-based parameterisa- 
tions of irreversible mixing processes for use in numerical 
ocean climate models. 

At a fundamental level, turbulent molecular diffusion in 
stratified fluids is important for at least two distinct — al- 
though inter-related — reasons: 1) for transporting heat dif- 
fusively across isopycnal surfaces - a process often referred 
to as 'diapycnal mixing'; 2) for dissipating available poten- 
tial energy, which contributes for a significant fraction — of- 
ten called the mixing efficiency — of the total dissipation 
of available mechanical energy ME, i.e., the sum of total ki- 
netic energy KE and available potential energy APE, which 
are defined by: 



KE^ j^P^dV, 



APE= I p{gz + I)dV 
'v 



p{gzr + Ir)dV, 



(1) 



(2) 



PE 



PEr 



where p is the density, v ~ {u, v, w) is the three-dimensional 
velocity vector, g is the acceleration of gravity, z is the verti- 
cal coordinate increasing upwar d, and / the spe cific internal 
energy. The APE is defined as in lLorenzl(ll955h as the differ- 
ence between the potential energy PE of the fluid (i.e., the 
sum of the gravitational potential energy GPE plus internal 
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energy IE) minus the potential energy PEr of a reference 
state that is the state of minimum potential energy achiev- 
able in an adiabatic re-ar rangement of the fluid parcels. As 
shown by IWinters & all fl995.) . the APE and PEr play a 
fundamental in the modern theory of turbulent mixing ow- 
ing to the fact that by construction PE,. is only affected by 
irreversible processes; as a result, measuring the time evolu- 
tion of the reference state provides a direct and objective way 
to quantify the amount of irreversible mixing taking place 
during turbulent mixing events, which is now commonly ex- 
ploited to diagnose mixin g in numerical experiments, e.g., 
Peltier & Caulfield I (12003V 

In the oceans, turbulent diapycnal mixing is required to 
transfer heat downward from the surface at a sufficiently 
rapid rate to balance the cooling of the deep ocean by high- 
latitudes deep water formation. In the oceanographic litera- 
ture, the most widely used approach to parameterise the ver- 
tical (diapy cnal) eddy diffusiv i ty Kn is based on the Osborn- 
Cox model (lOsborn and Cox I (119721) ): 



'^mixing 



N2 



iV2 



(3) 



which expresses Kp in terms of either the turbulent viscous 
kinetic energy dissipation ek or turbulent diffusive dissi- 
pation of available potential energy ep, where N'^ is the 
squared buoyancy frequency, and ^mixing = £p /^k is the 
ratio of the APE to KE dis sipation, which is often ca l led th e 
'mixing efficiency', e.g. iLindborg & Brethouwer ] (l2008h . 



Expressin g Kn in terms of Ek app ears to have been fi rst pro- 
posed by ^ Lilly & al. I (119741) and IWeinstock I ('1978') in the 
context of stratospheric turbulent mixin g, and adapted to the 



oceanographic case by lOsbornI (Il980h . The definition of 



mixing efficiency as a dissipations ratio adopted i n this pa- 
per appears to have been first proposed bv lOakevI (Il982h . 

Since both e p and ek are linked to the dissipation of me- 
chanical energy of which KE and APE represent the two 
main dynamically important forms, Eq. ^ makes it clear 
that turbulent diapycnal mixing is directly related to the me- 
chanical energy input in the oceans, but this link has been so 
far very rarely exploited in numerical ocean models. Rather, 
Kp is often regarded as a tunable parameter whose value 
is adjusted to reproduce the main observed features of the 



oceani c stratification. Such an approach was used bv lMunk 
d 19661) . who assumed the stratification to obey the vertical 
advective/diffusive balance: 



de d 



dz 



dz \ ^ dz 



(4) 



where is the potential temperature, and w the vertical ve- 
locity. Physically, Eq. (|4]i states that the upward advec- 
tion of cold water is balanced by the downward turbulent 
diffusion of heat, the rate of upwelling being set up by 
the rate of deep water formation. By using Eg. dill as a 



model for stratification profiles in the Pacific, iMunk I (119661) 



concluded that the canonical value K, 



10- 



Ys was 



apparently needed to explain the observed structure of the 
oceanic therm ocline. Subsequently, however, the validity of 
Munkl (ll966l) 's approach was questioned, as several observa- 
tional studies found Kp in the ocean interior to be typically 
sm aller by an order of m agnitude than Mu nk's value, e.g ., 
see ' Ledwell & all (11998 *) and the review by lOreggl (Il987l) . 
However, it seems widely re cognised today that K^ is highly 
variable spatially, prompting lMunk & Wunsch I (Il998h to re- 
interpret the value Kp = 10~^m^.s~^ as resulting from the 
overall effect of weak interior values combined with intense 
turbulent mixing in coastal areas or over rough topography. 

While the above approach is useful, it does not exploit 
the link between Kp and the mechanical sources of stir- 
ring sugg ested by Eq. dH. Clari fying this link was pio- 
neered by iMunk & Wunsch I dl998h . who translated the ad- 
vection/diffusion balance into one for the gravitational po- 
tential energy budget, which they argue must be a balance 
between the rate of GPE loss due to cooling and the rate of 
GPE increase due to turbulent diffusive mixing, i.e.. 



(5) 



d{GPEr) 




d{GPEr) 




dt 


cooling 


dt 


mixing 



this result being obtained by multiplying Eq. dUl by agpogz, 
after some manipulation involving integration by parts and 
the neglect of surface heating, where ae is the thermal expan- 
sion, g the acceleration of gravity, po a reference density, and 
z the vertical coordinate pointing upward. The subscript r 



is added here because it can be shown that Munk & Wunsch 



( Il998!) must actually pertain to the background G PEr bud 
Bjet, r ather than the AGPE budget, as shown by iTailleux 



get, n 
d2009 ^ 



d2009). This follows from the fact that cooling and turbulent 
molecular diffusion act as a GPE sink and source only for 
the background GPEr, as it is the opposite that holds for 
AGPE. If one assumes that density is primarily controlled 
by temperature for simplicity, the effect of mixing on GPEr 
is thus given by: 



d{GPEr) 



mixing 



dt 



poKpN' 1 



d ( d9 ^ 
Poctegz— I K„— I dV 



' dV, 
oz I 



(6) 



by using the result that N"^ = aggdO/dz in absence of 
salinity effects, and by assuming z = at the ocean sur- 
face, and no flux t hrough the ocean bottom. In their paper, 
Munk & Wunsch I dl998 i) neglected the nonlinearities of the 
equation of state, which amounts to regard ag as constant, in 
which case the above expression becomes: 



d{GPEr) 






dt 


mixing 


L 



PoKpN^dV. 



(7) 
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By using Eq. assuming ^mixing constant, this formula 
can be rewritten as follows: 



PoSpdV = 7„ 



d{GPEr) 


-1 


dt 


mixing 



D(APE) 




where D{APE) and D{KE) are the total volume-integrated 
diffusive dissipation of available potential energy and vis- 
cous dissipation of kinet ic energy respectively. To conclude, 
Munk & Wunschl(ll998l) linked the dissipation to production 
terms by assuming the balance D{KE) = G{KE), where 
G{KE) is the work rate done by the mechanical forcing due 
to the winds and tides. As a result, the above formula yields: 



G{KE) 



1 



'^mixing 



d{GPE) 



dt 



(9) 



cooling 



By estimating the rate of GPE loss due to cooling to 
be .4 TW, and by us i ng the canonical value ^mixing = 
0.2, lMunk & WunschI (Il998h concluded that G{KE) = 
0(2TW) of mechanical energy input was required to sus- 
tain turbulent diapycnal mixing in the oceans. Since the work 
of the wind stress against the sur face geostrophic velocity 
is widely agreed to be 0(1 TW), ' Munk & Wunsch I d 19981) 
suggested that the shortfall should be explained by the work 
rate done by the tides. The issue remains controversial, how- 
ever, because the role of the surface buoyancy forcing is not 
sufficiently well understood, as discussed in Tailleux ( 2009) . 

Another important issue in assessing the uncertainties 
associated with Eq. (|9]l concerns the importance of 
the nonlinearities of th e equation of state, neglected by 
Munk & Wunsch I (Il998h . As is well known, the nonlinear- 
ities of the equation of state are mostly responsible for the 
fluid "contracting upon mixing". This contraction is respon- 
sible for the actual increase in GPE^ due to mixing to be 
less than for a linear equation of state. In Eq. (|6]l, this can 
be seen from the fact that dag jdz is usually positive for a 
stably stratified fluid. Since z is negative by assumption, it 
follows that a correction factor is required that modifies Eq. 
(HJl as follows: 



d{GPE) 



dt 



= (1 - C,n)D{APE) 



mixing 

= (1 - CmhmixingDiKE). 

where Cm > 0, which in turn modifies Munk & Wunsch 



(10) 



(Il998l) 's constraint (Eq. ^) as foflows: 



G{KE) 



1 



^7: 



mixing 



d{GPE) 



dt 



(11) 



cooling 



where £ = 1 — Cm < 1- Based on the above arguments, 
Munk & Wunsch ( 1998b 's results are expected to underesti- 



point was first pointed out by iGnanadesikan & al. I (120051) 
who emphasised the importance of cabelling. Discussing the 
value of Cm or ^ to be used in Eq. ([9]) is beyond the scope 
of this paper. Note, however, that it is possible to construct 
stratifications with ^ not only smaller than one but al so pos - 



sibly even negative, as discussed by F ofonoff I (1199811200 Ih . 
The latter cases are interesting, because they are such that 
GPEr decreases upon mixing, not increases, in contrast to 
what is usually assumed. 

The main reason that GPEr is often assumed to increase 
as a result of turbulent mixing stems from that for an incom- 
pressible fluid with a linear equation of state, Eq. (O states 
that: 



d{GPE) 



dt 



D{APE), 



(12) 



i.e., that GPEr increases at the same rate that APE de- 
creases, which is classically interpreted as implying that the 
diffusively dissipa ted APE must be irr ev ersibly converted 



into GPEr, e.g., IWinters & al.1 (Il995h . iTailleuxl (l2009l) 



pointed out, however, that Eq. (fTzl i is at best only a good ap- 
proximation, not a true equality, since in reality the rates of 
GPEr increase and APE decrease are never exactly equal, 
and sometimes even widely different, because of the nonUn- 
ear character of the equation of state. 

In order to better understand how the net change in GPEr 
correlates with the total amount of APE diffusively dissi- 
pated during an irreversible turbulent mixing event, it is use- 
ful to examine the process of turbulent mixing in the light of 
classical thermodynamic transformations. To make progress, 
the conditions under which the diffusive exchange of heat be- 
tween fluid parcels takes place need to be known, but in prac- 
tice this is problematic, for it would require solving the full 
compressible non-hydrostatic Navier-Stokes equations down 
to the diffusive scales. Fortunately, it is often the case that 
stratified fluids at low Mach numbers are close to hydrostatic 
equilibrium, suggesting that the diffusive heat exchange be- 
tween parcels may reasonably be assumed to occur at ap- 
proximately constant pressure. If so, irreversible diffusive 
mixing must then be close to be a process conserving the to- 
tal potential energy PE — APE+PEr of the system, which 
implies that any amount AAPE^i < of diffusively dissi- 
pated APE must be irreversibly converted into background 
PEr, viz.. 



APEr = -AAPEdiff > 



(13) 



The implications for the net change in GPEr can be deter- 
mined from the definitions PEr = GPEr + lEr, APE = 
AGPE + AIE, and IE ^ AIE + lEr, which imply: 



and 



mate the constraint on G{KE), to the extent that 7„i 

the rate of GPEr loss due to coohng can be kept fixed. This AG PEr = — A APE, 



diff 



AIEr 



(14) 
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where AIEr is the net change in background internal energy 
taking place during the irreversible mixing event. As a result, 
the quantities ^ and Cm previously defined become: 



AGPEr 

wape:^\ 



AIEr 



\^APEd^ff\ 



Cy-fl, 



AIEr 



\AAPE, 



diffl 



(15) 



(16) 



Eqs. ( fTSl l and (fT6] l are important, because they establish that 
the nonlinearities of the equation of state — which are re- 
sponsible for the temperature and pressure dependence of 
a — can give rise to internal energy changes AIEr com- 
parable in magnitude with AAPEdif f and AGPEr dur- 
ing a turbulent mixing event. Such large lEr changes must 
in turn be associated with potentially large compressibility 
effects whose work against the pressure field may al so ex - 
pected to be large, as first demonstrated bv lTailleuxIdiooi) . 
In other words, the above formula suggest that the nonlin- 
earities of the equation of state may give rise to significant 
non-Boussinesq effects. So far, however, most numerical 
ocean models still make the incompressible and Boussinesq 
approximations, while at the same time using some version 
of the nonlinear equation of state for seawater Such an ap- 
proach yields values of ^ and Cm that are predicted by Eq. 
(|6]l, but since those values ultimately derive from initially 
making the Boussinesq approximation, it is unclear whether 
they can take into account the nonlinear character of the 
equation of state in a fully consistent manner. 

In fact, even when the net change AIEr appears to be 
small or ne gligible, se emingly justifying the incompressible 
assumption. i Tailleuxl (2009) argues that compressible effects 
may still be large, because one may show that AIEr can be 
decomposed as follows: 



AIEr = AIE, 



exergy 



AIEn 



(17) 



where IEq and lE^xergy = lEr — IEq are two subcom- 
ponents of lEr called the 'dead' and 'exergy' components. 
Physically, IEq represents the internal energy of a notional 
thermodynamic equilibrium state of uniform temperature Tq, 



whereas lE^ 



represents the internal energy associated 



with the vertical stratificat i on of the reference state. An im- 
portant result of iTailleux" (12009 ') is that the net changes in 



and IEq are related at leading order to AGPEr 



IE, 



and AAPEdi //as follows: 



AGPEr 



-AIE, 



exergy ; 



AIEo ~ -AAPE, 



dtff 



> 0, 



(18) 

(19) 



to a very good approximation in a nearly incompressible fluid 
such as water or seawater. These relations state that turbulent 
molecular diffusion primarily dissipates APE into 'dead' in- 
ternal energy IEq, while simultaneously causing a transfer 



results in an increase of the equivalent thermodynamic tem- 
perature Tq, whereas the latter effect results in the smoothing 
out of dTr/ dz. This contrasts with the standard interpretation 
that turbulent molecular diff usion irreversib l y conv erts APE 
into GPEr, as proposed by IWinters & al. I (Il995b . The dif- 
ferences between the two interpretations are schematically 
illustrated in Fig. 1. The main reason why compressibility 
effects may be important even if AIEr ~ is because vol- 
ume changes are primarily determined by AIE exergy, not by 
AIEo or AIEr. 

As regards to the empirical determination of the mixing 
efficiency ^mixing — £p/£K, the above remarks are impor- 
tant because AGPEr and | AAPEdi // 1 are currently widely 
thought to physically represent the same quantity, prompting 
many studies to actually estimate e p from measuring the net 
changes in GPEr, e.g., JvlcEwaiL (1983a,b); B arry. (20011) . 
For the reasons discussed above, however, this makes sense 
only if ^ can be ascertained to be close to unity, as if not, 
the relevant value of ^ is then required. One of the main ob- 
jective of this paper is to establish that the behaviour of ^ is 
closely connected to the sign and amplitude of the following 
parameter: 



(20) 



between GPEr and IE, 



exergy 



. Physically, the former effect 



d I aP' 
dz \pCp^ 

where a is the thermal expansion coefficient, P is the pres- 
sure, p is the density, and Cp is the specific heat capac- 
ity at constant pressure, while salinity is assumed to be 
uniform throughout the domain. Physically, the parameter 
T = aP/{pCp) = PT/T, where T is the adiabatic lapse 
rate, represents the fraction of the amount of heat SQ re- 
ceived by a parcel in an isobaric process that can be converted 
into work. As a result, T is expected to be the main param- 
eter controlling the net change in GPEr due to the turbulent 
diffusive heat exchange between fluid parcels. 

From the viewpoint of turbulent mixing, the main diffi- 
culty posed by a nonlinear equation of state is to make it 
possible for different vertical stratification to share the same 
profile N{z) without necessarily having the same T{z) verti- 
cal profile. From a dynamical viewpoint, this is not expected 
to be a problem as long as the dynamical evolution of KE 
and APE, as wefl as D{KE) and D{APE), remain mostly 
controlled by N{z) at leading order, as is usually assumed. 
If so, the dissipations ratio D{APE)/ D{KE), and hence 
the bulk mixing efficiency ^mixing, can then be assumed to 
be unaffected by the nonlinearities of the equation of state at 
leading order. The main objective of this paper is to verify 
that D{APE) appears indeed to be largely insensitive to the 
T(z) vertical profile, and hence mostly controlled by N{z). 
If so, we can safely conclude that it must also be the case for 
D{KE), since there is even less reasons to believe that the 
latter could be affected by T{z). This could be directly ver- 
ified through direct numerical simulations of turbulent strat- 
ified mixing using a fully compressible Navier-Stokes equa- 
tions solver, which we hope to report on in the future. On 
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the other hand, the net change in GPEr is expected to be 
extremely sensitive to T(z). Most of the paper is devoted 
to verify that this is indeed the case, and to find ways to re- 
late the net change in GPEr to the sign and magnitude of 
dT/dz. Section 2 provides a theoretical formulation of the 
issue discussed. Section 3 discusses the methodology, while 
the results are presented in Section 4. Finally, section 5 sum- 
marises and discusses the results. 



2 Theoretical formulation of the problem 

2. 1 Energetics of mixing 

A key issue in the study of turbulent mixing is understanding 
the links betw een stirring and mixing. As first discussed by 
Eckart the two processes can be rigourously sepa- 

rated if one notes that the probability density functions (pdf 
in short) of the adiabatically conserved quantities (i.e., en- 
tropy and salt for seawater) are only affected by the irre- 
versible mixing due to the molecular diffusion of heat and 
salt, but not by the adiabatic shuffli ng of the parce ls due to 



the stirring process. The link with iLorenz I (Il955 l)'s avail 



able potential energy framework comes from the fact that 
Lorenz's reference state, i.e., the state whose potential en- 
ergy is minimised by an adiabatic re-arrangement of the fluid 
parcels, coincides with the above-mentioned pdf, which was 
exploited by iWinters & al. I (1 19951) to provide a new way to 
rigourously quantify irreversible mixing simply from diag- 
nosing the temporal evolution of the reference state. In a 
Boussinesq fluid with a linear equation of state, the role of 
entropy is played by either temperature or density. In the fol- 
lowing, the fluid will be assumed to be either freshwater or 
seawater with u niform salinity t hroug hout the fluid. 

As shown bv I Winters & al. I d 1995 1) (using somewhat dif- 
ferent notations), the energetics of freely decaying turbulence 
in an insulated domain is based on the following evolution 
equations for the volume-integrated kinetic energy (KE), 
available potential energy {APE), and background gravita- 
tional potential energy (GPEr): 



d{KE) 
dt 

d{APE) 
dt 

djGPEr] 

dt 



= ~C{KE, APE) - D{KE), 



C{KE,APE) - D{APE), 



Wr 



W, 



rAamtnar 



w, 



r.turhulent 



(21) 



(22) 



,(23) 



where C{KE, APE) is the so-called buoyancy flux, which 
physically represents the reversible conversion between KE 
and APE, while all other terms represent irreversible pro- 
cesses, with D{KE) denoting the viscous dissipation of 
KE, D{APE) the diffusive dissipation of APE, and 
the rate of change of GPEr due to molecular 



and laminar contribution. Note that the above equations are 
domain-averaged, not local formulations, which are expected 
to be well suited for understanding laboratory experiments of 
turbulent mixing for which lateral fluxes of APE and KE 
can be ignored. 

As discussed bv iTailleux I (12009'). Eqs. (EBEll) provide 
a unifying way to describe the energetics of both the incom- 
pressible Boussinesq and compressible Navier-Stokes equa- 
tions, by adapting the definitions of the energy reservoirs 
and energy conversion terms to the particular set of equa- 
tions considered. Expli cit expres s ions f or D{APE) and 
Wr^mixing SIC given by iTaiUeux I (l2009l) in the particular 
cases of: 1) a Boussinesq fluid with a linear and nonlin- 
ear equation of state in temperature; 2) for a compressible 
thermally-stratified fluid obeying the Navier-Stokes equa- 
tions of state with a general equation of state depending on 
temperature and pressure. These expressions are recalled 
further below for case 2). While Wr.iaminar is well under- 
stood to be a conversion between IE and GPEr, the na- 
ture of the energy conversions associated with D{APE) and 
Wr,turbuient IS Still a matter of debate. Currently, it is widely 
assumed that D{APE) and Wr,turbuient represent the same 
kind of energy conversion, namely the irreversible conver- 
sion of APE into GPEr owing to the fact that for a Boussi- 
nesq fluid with a linear equation of state (referred to as the 
L-Boussinesq model hereafter), one has the exac t equality 



D(AP E) = Wr,turbuient- It was pointed out by Tailleux 



(l2009h that this equality is a serendipitous artifact of the L- 
Boussinesq model, which does not hold for more accurate 
forms of the equations of motion. More generally, 'Tailleux I 
(2009) found that the ratio ^ Wr,turbuientl D{APE) is 
not only systematically lower than unity for water or sea- 
water, but can in fact also takes on negative values, as pre- 
viously discussed by .Fofonoff (1962, 1998, 2001 ) in a se- 
ries of little known papers. In other words, the equality 
D{APE) — Wr,turbuient IS Only a mathematical equality, 
not a physical equality, by defining a physical equality as a 
mathematical equality between two quantities that persists 
for the most accurate forms of the governing equations of 
motion. To clarify the issue, .Taille ux (2009) sought to un- 
derstand the links between D{APE), Wr,mixing and internal 
energy, by establishing the following equations: 



djlEg) 
dt 

d{IE, 



D{KE) + D{APE), 



exergy ) 



dt 



(24) 



(25) 



r^7mxtng 



diffusion, which is customarily decomposed into a turbulent 



which demonstrate that the viscously dissipated KE and dif- 
fusively dissipated APE both end up into the dead part of 
internal energy lEo, whereas Wr, mixing represent the con- 
version rate between GPEr and the 'exergy' component of 
internal energy lEexergy ■ A schematic energy flowchart il- 
lustrating the above points is provided in Fig. 1 . 
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2.2 Efficiency of mixing and mixing efficiency 



AGPEr = Wr 



Wr.t 



urhulent 



r .Laminar 



, (29) 



The APE framework introduced bv lWinters & al. I ( Il995h for 
a Boussi nesq fluid with a linear equation of state, and ex- 



tended by lTailleux I (l2009h to a fully compressible thermally- 



stratified fluid, greatly simplifies the theoretical discussion of 
the concept of mixing efficiency. To that end, it is useful to 
start with the evolution equation for the total "available" me- 
chanical energy ME ~ KE + APE, obtained by summing 
the evolution equations for KE and APE, leading to: 



= - [D{KE)+D{APE)]. 



(26) 



d{ME) 
Jt 

Eq. (|26] |, along with Eq. (|24] |. are very important, for 
they show that both viscous and diffusive processes con- 
tribute to the dissipation of ME into deal internal energy 
IEq. From this viewpoint, understanding turbulent diapyc- 
nal mixing amounts to understanding what controls the ratio 
Iraixmg = D{APE) /D{KE), that is, the fraction of the to- 
tal available mechanical energy dissipated by molecular dif- 
fusion rather than by molecular viscosity. The amount of ME 
dissipated by molecular diffusion, i.e., D{APE), is impor- 
tant, because it is directly related to the definition of turbulent 
diapycnal diffusivity, as said above in relation with Eq. (O. 

The link between the dissipation mixing efficiency and 
more traditional definitions of mixing efficiency can be clari- 
fied in the light of the above energy equations, by investigat- 
ing the energy budget of a notional "turbulent mixing event", 
defined here as an episode of intense mixing followed and 
preceded by laminar conditions (i.e., characterised by very 
weak mixing), during which KE and APE undergo a net 
change change AKE < and AAPE < 0. As far as we 
understand the problem, most familiar definitions of mixing 
efficiency appear to implicitly assume AAPE « 0, as is the 
case for a turbulent mixing event developing from a unsta- 
ble stratified shear flow for instance, e.g., P eltier & Caulfield . 



(120031) . This point can be further clarified by comparing the 
energetics of turbulent mixing events developing from the 
shear flow instability with that developing from the Rayleigh- 
Taylor instabiUty, treated next, which by contrast can be re- 
garded as having the idealised signature AKE « and 
AAPE < 0. 

In the case of the stratified shear flow instability, assumed 
to be such that AKE < and AAPE « 0, integrating flie 
above energy equations over the time interval over which the 
turbulent mixing event takes placeQ yields: 



AKE = ~C{KE,APE) - D{KE), 



(27) 
(28) 



= C{KE, APE) - D{APE), 

'it is usually assumed that the time average should be short 
enough that the viscous dissipation of the mean flow can be ne- 
glected. Alternatively, one should try to separate the laminar from 
the turbulent viscous dissipation rate. The following derivations as- 
sume that the viscous dissipation is dominated by the dissipation of 
the turbulent kinetic energy rather than that of the mean flow. 



where the overbar denotes the time integral over the mix- 
ing event. For a Boussi nesq fluid with a linear equation 
of state, I Winters & al. I (1 19951) showed fliat D{APE) = 
W r.turbuient- If wc Combine the latter result with the APE 
budget (i.e., Eq. (|28]|). one sees that one has the triple equal- 
ity: 



C{KE,APE) ^ D{APE) = W 



r.turhulent ■ 



(30) 



The triple equality Eq. ( l30b suggests that any of the three 
quantities C{KE, APE), D{APE), or Wr,turbuient can a 
priori serve to measure "the fraction of the kinetic energy that 
appears as the potential energy of the stratification", which 
is the traditio nal definition of the flux Richardson number 
proposed by Linden (Il979l) . Historically, the buoyancy flux 



C{KE,APE) is the one that was initially regarded as the 
natural quantity to use for that purpose in an overwhelming 
majority of past studies of turbulent mixing. As a result, most 
existing studies of turbulent mixing define the turbulent di- 
apycnal diffusivity, mixing efficiency, and flux Richardson 
number in terms of the buoyancy flux as follows: 



^ C{KE,APE) 



flux 



C{KE,APE) 
D{KE) 



R 



flux 



C{KE, APE) 



C{KE, APE) + D{KE) 



(31) 



(32) 



(33) 



It is easily verified that the above equa tions are consistent 
with those considered bv lOsborn (Il980l) for instance. Phys- 
ically, however, there are fundamental problems in using the 
buoyancy flux to quant ify irreversible diffu s ive rn i xing, be- 
cause as po inted out by Caulfield & Peltier I (I2OO0I) . IStaquet 



(I2OOOI) and lPeltier & Caulfield I (l2003l) . CiKE. APE) rep 



resents a reversible energy conversion, which usually takes 
on both large positive and negative values before settling on 
its long term average D{APE). Moreover, as pointed out be- 
low, the buoyancy flux is only related to irreversible diffusive 
mixing only if AAPE w holds to a good approximation, 
for otherwise, it becomes also related to the irreversible vis- 
cous dissipation rate as shown by the KE budget (Eq. (l27Ti). 
Eq. dSOl l makes it possible, however, to use either D{APE) 
or W r,turbuient instead of C{KE, AP E) in the definitions 
(ISTTl a nd (l32b. For this reason, both ICaulfield & Peltier I 
( 2000l) and Staquet ( 2000l) proposed to measure the effi- 
ciency of mixing based on Wr,turbuient, i-e.. 



GPEr 



K 



GPEr 

1 mixing 



w 



r.turhulent 



Wr 



turbulent 



D{KE) 



(34) 



(35) 
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Wr.t 



urhulent 



Wr. 



turbulent 



D{KE) ' 



(36) 



such a definition being motivated bv lWinters & al. ( 1995 )'s 
interpretation that D{APE) and Wr ^turbulent represent the 
same energy conversion whereby the diffusively dissipated 
APE is irreversibly converted into GPEj.- The parame- 
was called the " cumulative mixing efficiency" 



ter R 



GPEr- 

■s- 



bv lPeitier & Caulfield ' ('2003') and modified flux Richardson 
n umber by|S taquet (2000). As argued in Tailleux (200i), it 
is D{APE), rather than W r, turbulent, that directly measures 
the amount of KE eventually dissipated by molecular dif- 
fusion via its conversion into APE, suggesting that the flux 
Richardson number should actually be defined as: 



jyDAPE 



D{APE) 



D{KE) + D{APE) 



(37) 



While the above formula makes it clear that all above defini- 
tions of Rf are equivalent in the particular case considered, 
it is easily realized that they will in general yield different 
numbers if one relaxes the assumption AAPE « in Eq. 
( |28] ), as well as the assumption of a linear equation of state, 
yielding a ratio ^ = WrMirbuient / D{APE) that is gener- 
ally lower than unity and sometimes even negative for water 
or seawater For this reason, it is crucial to understand the 
physics of mixing efficiency at the most fundamental level. 
From the literature, it seems clear that most investigators's 
idea about the flux Richardson number is as a quantity com- 
prised between and 1. From that viewpoint, the dissipation 
flux Richardson number R^^^^ is the only quantity that sat- 
isfies this property under the most general circumstances, as 
cases can easily be constructed for which both Wr,turbuient 
and C{KE,APE) are negative. Indeed, cases for which 
^ < are described in this paper, whereas C{KE, APE) is 
easily shown to be negative in the case of a turbulent mixing 
event for which all mechanical energy is initially provided 
entirely in APE form. In that case, assuming AAPE < 
and AKE « in the above energy budget equations yields: 



C{KE, APE) = D{APE) + AAPE = -D{KE), (38) 



which shows that this time, C{KE, APE) directly measures 
the amount of viscously dissipated kinetic energy, rather 
than diapycnal mixing. The latter case is relevant to under- 
sta nd the energy budget of the Rayleigh-Taylor instability, 
see lDalziel & al I (i2008) for a recent discussion of the latter 



2.3 Link between D{APE) and Wr^ 



7mxtng 



In order to help the reader understand or appreciate why the 
ratio £, = Wr,turbuient / D{APE) is generally lower than 
unity for water or seawater, and hence potentially signif- 
icantly different from the predictions of the L-Boussinesq 



model, it is useful to examine the structure of Wr.i 



and 



D{APE) in more details. As shown by Tailleux ( 20091) . 



the analytical formula for the latter quantities in a fully com- 
pressible thermally-stratified fluid are given by: 



D{APE) 



n P 

" '-V- (KpCpVT)dF, 



PrCpr 
T-Tr 



V • [npCpVT) dV, 



(39) 



(40) 



where as before a is the thermal expansion coefficient, P is 
the pressure, Cp is the specific heat capacity at constant pres- 
sure, p is density, with the subscript r indicating that values 
have to be estimated in their reference state. The parame- 
ter T — aP/{pCp) plays an important role in the problem. 
Physically, it can be shown that in an isobaric process during 
which the enthalpy of the fluid parcel increases by dH, the 
parameter T represents the fraction of dH that is not con- 
verted into internal energy, i.e., the fraction going into work 
(and hence contributing ultimately to the overall net change 
in GPEr). As a result, T plays the role of a Carnot-like ther- 
modynamic efficiency. In Eq. (l39l l. denotes the value that 
T would have if the corresponding fluid parcel was displaced 
adiabatically to its reference position. 

In order to compare these two quantities, we expand T as 
a Taylor series around P = Pr, viz.. 



T = Tr+Tr{P - Pr 



(41) 



where Tr — arTr/ [prCpr) is the adiabatic lapse rate. At 
leading order, therefore, one may rewrite D{APE) as fol- 
lows: 



D{APE) = / 
Jv 



ar{Pr - P) Tr 



prCj 



pr 



T 



V • {KpCpVT)dV + 



Wr, 



UrTrP 
Pr Cpr T 



{Tr - T) arPr 

(upCpVT) dV + • 



V • {npCpVT) dV 



(42) 



These formula shows that D{APE) can be written as the 
plus some corrective terms . One sees that 



sum of W, 



r^mixtng 



the L -Boussinesq model's results derived bv [Winters & al. 
(Il995i) can be recovered in the limit T « Tr, P « —pogz, 
ar/{prCpr) ~ Q!o/(/3oC'po), pCp ~ PoCpo, whcrc the 
subscript refers to a constant reference Boussinesq value, 
yielding: 



D{APE) « Wr,, 



Wr 



.amtnar 



Wr 



r. turbulent ■ 



(43) 



These results, therefore, demonstrate that the strong corre- 
lation between D{APE) and Wr, mixing originates in both 
terms depending on molecular diffusion in a related, but nev- 
ertheless distinct, way, the differences between the two quan- 
tities being minimal for a linear equation of state. The fact 
that the two t erms are never exactl y equal in a real fluid 
clearly refutes Iwinters & al. ( 1995 )'s widespread interpre- 
tation that D{APE) and Wr, turbulent physically represents 
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the same energy conversion whereby the diffusively dissi- 
pated APE is irreversibly converted into GPEr- In real- 
ity, D{APE) and Wr.turbuient represent two distinct types 
of energy conversions that happen to be both controlled by 
stirring and molecular diffusion in related ways, which ex- 
plains why they appear to be always strongly correlated, and 
even exactly equal in the idealised limit of the L-Boussinesq 
model. If one accepts the above point, then it should be clear 
that what is now required to make progress is the understand- 
ing of what controls the behaviour of the parameter ^, since 
the knowledge of the latter is obviously crucial to make infer- 
ences about turbulent diapycnal mixing from measuring the 
net changes of GPE^ for instance. The purpose of the nu- 
merical simulations described next is to help gaining insights 
into what controls ^. 



3 Methodology 

To get insights into how the equation of state of seawa- 
ter affects turbulent mixing, we compared D{APE) and 
Wr.turbuient for a uumbcr of different stratifications having 
the same buoyancy frequency vertical profile N, but different 
vertical profiles with regard to the parameter aP/{pCp), as 
illustrated in Fig. 2. The quantities D{APE) and Wr^mixing 
were estimated from Eqs. ( [39] l and (l40l i. while Wr,turbuient 
was estimated from 



.turbulent 



mixing 



r, laminar •, 



(44) 



where Wr.iaminar was obtained by taking T = Tr in 
the expression for Wr. mixing- The quantities D{APE) 
and Wr,turbuient wcrc estimated numerically for a two- 
dimensional square domain discretised equally in the hori- 
zontal and vertical direction. In total, 27 different stratifica- 
tions were considered, all possessing the same squared buoy- 
ancy frequency N"^ illustrated in the left panel of Figure 2, 
but different mean temperature, salinity, and pressure result- 
ing in different profiles for the aP/{pCp) parameter illus- 
trated in the right panel of Figure 2. In all cases considered, 
the pressure varied from Pmin to Pmax = Pmin + lOdhar, 
with Pmin taking the three values (0 dbar, 1000 dbar, 2000 
dbar). In all cases, the salinity was assumed to be con- 
stant, and taking one of the three possible values S = (30 
Psu, 35 psu, 40 psu). With regard to the temperature pro- 
file, it was determined by imposing the particular value 
Tmax = T{Pmin) at the top of the fluid, with all remain- 
ing values determined by inversion of the buoyancy fre- 
quency N"^ common to all profiles by an iterative method. 
The imposition of a fixed buoyancy profile N, salinity S, 
pressure range, and minimum temperature Tmin was found 
to yield widely different top-bottom temperature differences 
T {Pmin) ~T {Pmax), ranging from a few tenths of degrees to 
about 4 degrees C depending on the case considered, as seen 
in Fig. 3. In each case, the thermodynamic properties of 
the fluid were estimated from the Gibbs function of .Feistel ■ 



(120031) . Specific details for the temperature, pressure, and 
salinity in each of the 27 experiments can be found in Table 
1 along with other key quantities discussed below. 

Numerically, the two-dimensional domain used to quan- 
tify D{APE) and Wr, turbulent was discretised into Npi x 
Npj points in the horizontal and vertical, with Npi — Npj = 
100. Mass conserving coordinates were chosen in the verti- 
cal, and regular spatial Cartesian coordinate in the horizontal. 
For practical purposes, the vertical mass conserving coordi- 
nate can be regarded as standard height z, as the differences 
between the two types of coordinates were found to be in- 
significant in the present context, and thus chose Ax = Az. 
In order to compute D{APE) and Wr.turbuient for turbulent 
conditions, we modelled the stirring process by randomly 
shuffling the fluid parcels adiabatically from resting initial 
conditions. Shuffling the parcels in such a way requires a 
certain amount of stirring energy, which is equal to the avail- 
able potential energy APE of the randomly shuffled state. 



4 Results 

For each of the 27 particular reference stratifications consid- 
ered, synthetic turbulent states were constructed by gener- 
ating hundreds of random permutations of the fluid parcels, 
thus simulating the effect of adiabatic shuffling by the stirring 
process, in each case yielding a particular value of D{APE), 
Wr, mixing, Wr, turbulent and APE. One Way to illustrate 
that Wr,turbuient depends more sensitively on the equation 
of state than D{APE) is by plotting each quantity as a func- 
tion of APE, as illustrated in Fig. 4. Interestingly, the figure 
shows that all values of D{APE) appear to be close to a lin- 
ear straight line, with no obvious sensitivity to the particular 
value of T. In contrast, the right panel of Fig. 4 demonstrates 
the sensitivity of Wr,turbuient to T, as a separate curve is ob- 
tained for each different stratification. Note that one should 
not construe from Fig. 4 that D{APE) is a linear function 
of APE. Physically, D{APE) depends both on the APE, 
as well as on the spectrum of the temperature field. It so 
happens that the method used to randomly shuffle the parcels 
tends to artificially concentrate all the power spectrum at the 
highest wavenumbers, the effect of which being to suppress 
one degree of freedom to the problem, which is responsible 
for the appearance of a linear relationship between D {APE) 
and APE in Fig. 4. It is easy to convince oneself, however, 
that stratifications can be constructed which have the same 
value of APE, but widely different values of D{APE). 

In order to understand how the equation of state affects 
Wr,turbuient, it is uscful to rewrite Wr,niixing as givcu by 
Eq. ( [39] ) as foUows: 



Wr,, 



npCpVT ■ V 

Y \Pr^pr 



dV 



f d f arPr \ dTr„ „n 

- ^"^^97L^hr"^"^"'^^+'" ^^^^ 

Jy OZr xPr'^pr J 'J'^r 
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by using an integration by parts, assuming insulated bound- 
aries, and using the approximation VT w VT^ + 0(T — Tr), 
by noting that the reference quantities depend only upon z^. 
Eq. gill suggests that Wr,mtxing and Wr,turbuient are pri- 
marily controlled by the vertical gradient of T = aP/{pCp), 
and that both Wr^^xing and Wr^t^^rbuient are likely to be pos- 
itive only when dT /dz is negative. This is obviously the case 
when the vertical variations of a/ (pCp) can be neglected, as 
in this case dT /dz « a/ {pCp)dP / dz « —ag/Cp < 0, as- 
suming the pressure to be hydrostatic. The case when the 
vertical gra dient of aP/jpCr,) is positiv e was extensively 
discussed by iFofonoff I i 1 962l 1 1 99811200 11) . and can be easily 
encountered in the oceans. 

In all experiments considered, we found the ratio ^ = 
Wr,turbuientl D{APE) to be Systematic ally l ower than 
unity, as already pointed out in TaiUeuxl (l2009l) . In order 
to better understand how dT /dz controls the behaviour of 
Wr,turbuient, the ratio ^ = D{APE) /Wr^turbuient was av- 
eraged over all randomly shuffled states separately for each 
stratification, the results being summarised in Fig. 5 and Ta- 
ble 1, along with the minimum value of dT /dz, as well as 
with the top-bottom difference AT = T{Pmin) — ^{Pmax)- 
Panels (a) and (c) show that as long that dT < 0, the equality 
Wr,turbuient ~ D{APE) holds to a rather good approxima- 
tion, up to a factor of 2, the approximation being degraded 
at the lowest temperature and salinity. Note, however, that 
in the cases considered, ^ > only at atmospheric pressure, 
with ^ being systematically negative at Pmm — lOOOdbars 
and Pmin — 2000 dbar respectively. Both Table 1 and 
Fig. 5 (a) and (c) show that ^ becomes increasingly neg- 
ative as [dT/dz]min becomes increasingly large and posi- 
tive, the worst case being achieved for the lowest T, lowest 
salinity, and highest pressure. As a further attempt to un- 
derstand this behaviour, we also computed the average ratio 
AGP E / APE for each particular reference stratification. In- 
terestingly, we find that the classical case ^ « 1 coincide with 
APE K, AGPE, as expected in the Boussinesq approxima- 
tion. We find, however, that the decrease in ^ coincides with 
AGPE being an increasingly bad approximation of APE. 
As the latter implies that AIE becomes increasingly impor- 
tant, it also implies that compressible effects become increas- 
ingly important. This suggests, therefore, that the effects of a 
nonlinear equation of state are apparently strongly connected 
to non-Boussinesq effects, a topic for future exploration. 

The key point of the present results is that while there 
exist stratifications such that Wr^turbuient ~ D{APE) to 
a good approximation, and hence that conform to classical 
ideas about turbulent mixing in a Boussinesq fluid with a lin- 
ear equation of state, there also exist stratification for which 
Wr.turbuient and D{APE) differ radically from each other. 
The main reason why this is not more widely appreciated 
is suggested by the results summarised in Table 1, which 
shows that Wr, turbulent ~ D{APE) appears to hold well 
under normal temperature and pressure conditions, which are 
usually those encountered in most laboratory experiments 



of turbulent mixing. In that case, the classical results of 
Boussinesq theory are applicable, and there is no problems in 
measuring the mixing efficiency of turbulent mixing events 
from measuri ng the net change in GPE^, as often done, e.g., 
Barry I (1200 ih . in acc ordance with the definition of m ixing ef- 
ficiency proposed by Caulfield & Peltier ( 2000l) and Staquetl 
(.2000.) . since ^ « 1 to a good approximation. Temperature, 
salinity, and pressure conditions in the real oceans can be 
very different than in the laboratory, however, especially in 
the abyss. In the latter case, the present results suggest not 
only that £, can potentially become very large and negative, 
but that the discrepancy between AGPE and APE can be- 
come significant to the point of making the Boussinesq ap- 
proximation and the neglect of compressible effects very in- 
accurate. This point seems important in view of the current 
intense research effort devoted to understanding tidal mixing 
in the abyssal ocean s that was prompte d a dec ade ago by the 
influential study by Munk & Wunsch I (1998). The point is 
also important because values of mixing efficiency published 
in the literature have been traditionally been reported without 
mentioning the associated value of ^, which may explain part 
of the spread in the published values, and adds to the uncer- 
tainty surrounding this crucial parameter The present results 
suggest that an important project would be to seek to recon- 
struct the missing values of ^, which is in principle possible if 
sufficient information about the ambient conditions are avail- 
able. 



5 Conclusions 

The nonlinearities of the equation of state for water or sea- 
water make it possible for a stratification with given mean 
vertical buoyancy profile N to have widely different ver- 
tical profiles of the parameter T — aP/{pCp), depend- 
ing on particular oceanic circumstances. The main result 
of this paper is that the sign and magnitude of dT/dz 
greatly affect Wr^turbuient — the turbulent rate of change of 
GPEr — while they correspondingly little affect D{APE), 
the dissipation rate of APE. As a result, the ratio ^ = 
Wr^turbuient / D{APE) is in general lower than unity, and 
sometimes even negative, for water or seawater. For this rea- 
son, the fact that D{APE) and Wr,turbuient happen to be 
identical for a Boussinesq fluid with a linear equation of state 
appears to be a very special case, which is rather misleading 
in that it fails to correctly address the wide range of values 
assumed by the parameter ^ in the actual oceans, while also 
leading to the widespread erroneous idea that the diffusively 
dissipated APE is irreversibly converted into GPEr, and 
hence that turbulent mixing always increase GPE. As far as 



we understand the problem, based on the analysis of lTailleux 
(2009), D{APE) and Wr,turbuient represent two physically 
distinct kinds of energy conversion, the former associated 
with the dissipation of APE into 'dead' internal energy, and 
the latter associated with the conversion between GPEr and 
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the 'exergy' part of internal energy. The former is always 
positive, while the latter can take on both signs, depending 
on the particular stratification. 

From the viewpoint of turbulence theory, the present re- 
sults indicate that the equality D{APE) = Wr^turbuient 
obtained in the con text of the L-Boussinesq model by 



Winters & al. I (Il995h should only be construed as implying 
a strong correlation between D{APE) and Wr^turbuient, not 
as an indication that the diffusively dissipated APE is con- 
verted into GPEr- As the present results show, the correla- 
tion between the two rates strongly depends on the nonlin- 
earities of the equation of state. Fundamentally, D{APE) 
and Wr,turbuient appear to be correlated because they both 
depend on molecular diffusion, and on the gradient of the 
adiabatic displacement C, — z — Zj. of the isothermal surfaces 
from their reference positions. Based on the present results, 
the ratio ^ = Wr.turbuient/ D{APE) appears to be deter- 
mined at leading order mostly by the sign and magnitude of 
dT/dz ~ d/dz[aP/ {pCp). Further work is required, how- 
ever, to clarify the precise link between ^ and dT /dz under 
the most general circumstances, which will be reported in a 
subsequent paper. 

The present results are important, because they show that 
the two following ways of defining a flux Richardson number 
Rf and mixing efficiency '^mixing, viz.. 



..DAPE 



D{APE) 
D{KE) '■ 



r,j,pE ^ DjAPE) 

i D{APE) + D{KE) 



(46) 



(47) 



called the dissipation mixin g efficiency and flux Richardson 
number by Tailleux ( 2QQ9b . and 



GPEr ^^r, turbulent 
Imixing 



D{KE) 



R 



GPEr 
f 



Wr. 



turbulent 



r.turbulent 



+ D{KE) ' 



(48) 



(49) 



as proposed by ICaulfield & Peltier ] (l2000h and IStaquet 
lloGO), which are equivalent in the context of the L- 
Boussinesq model, happen to be different in the context of 
a real compressible fluid, as the conversion rules 



GPEr ^f^DAPE 
Imixing S Imixing "! 



R 



GPEr 



^R 



DAPE 



1 - (1 - ■ 



(50) 



(51) 



now involve the parameter ^. Note that historically the flux 
Richardson number was defined bv lLindenI (1 19791) as "The 



fraction of the kinetic energy which appears as the poten- 
tial energy of the stratification." Physically, the kinetic en- 
ergy that appears as the potential energy of the stratification 
is the fraction of kinetic energy being converted into APE 



and ultimately dissipated by molecular diffusion. This frac- 
tion is therefore measured by D{APE), not by Wr^turbuient, 
since the latter technically represents the "mechanically- 
controlled" frac tion of int e rnal e nergy converted into GPEr, 
if one accepts iTaiUeux I (l2009h 's conclusions. From this 



viewpoint, it is R? ^^^ rather tha n i??^^'' that appears to 



I 1 J ^ 

be consistent with iLinden I (Il979h 's definition of the flux 



Richardson number, a nd hence 2£^nq rather than 



..GPE^ 
mixing 



that is consistent with lOsborn I (Il980l) 's definition of mixing 
efficiency. 

From a practical viewpoint, however, the above conceptual 
objections against 7,^ ■ ^ ■ ""^ and Rf do not mean that it is 
equally physically objectionable to seek estimating the effi- 
ciency of mixing from measuring the net changes in GPEr 
taking plac e during a turb ulent mixing event, as is commonly 
done, e.g.. lBarrvl (1200 li) . Such a method is perfectly valid, 
owing to the correlation between D{APE) and Wr,turbuient ■ 
The present results show, however, that such an approach re- 
quires the knowledge of the parameter ^, which is usually 
not supplied. For most laboratory experiments performed at 
atmospheric pressure, the issue is probably unimportant, as ^ 
appears to be generally close to unity in that case. The issue 
becomes more problematic, however, for measurements car- 
ried out in the ocean interior, as there is less reason to assume 
that ^ ~ 1 will be necessarily verified. A critical review of 
published values of jmixing would be of interest, in order to 
identify the cases potentially affected by a value of ^ signifi- 
cantly different from unity. 

So far, we have only considered the case of an equation of 
state depending on temperature and pressure only, by holding 
salinity constant. In practice, however, many studies of tur- 
bulent mixing are based on the use of compositionally strat- 
ified fluids. Understanding whether can be significantly 
different from unity in that case remains a topic for future 
study. 

Acknowledgements. The author acknowledges funding from the 
RAPID programme. Rainer Feistel, Michael Mclntyre and an 
anonymous referee are also thanked for provided comments that led 
to significant improvements in presentation and clarity. 



References 

Barry, M.E., Ivey, G. N., Winters, K. B., and Imberger, I. 
2001 Measurements of diapycnal diffusivities in stratified fluids. 
/ Fluid Mech. 442, 267-291. 

Caulfield, CP. c& Peltier, W.R. 2000 The anatomy of the mix- 
ing transition in homogeneous and stratified free shear layers. / 
Fluid Meek, 413, 1-47. 

Dalziel, S. B., Patterson, M. D., Caulfield, C. P., & 
COOMARASWAMY, LA. 2008 Mixing efficiency in high-aspect- 
ratio Rayleigh-Taylor experiments. Phys. Fluids, 20, 065106. 

ECKART, C. 1948 An analysis of the stirring and mixing processes 
in incompressible fluids. / Mar. Res., 1, 265-275. 

Feistel, R. 2003 A new extended Gibbs thermodynamic potential 
of seawater. Prog. Oceanogr. 58, 43-1 14. 



R. Tailleux: Is mixing efficiency affected by nonlinear equation of state? 



11 



FOFONOFF, N.P. 1962 Physical properties of seawaten M.N. Hill, t 

Ed., The Sea, Vol. 1, Wiley-Interscience, 3-30. 
FOFONOFF, N.P. 1998 Nonlinear limits to ocean thermal structure. 

J. Mar. Res. 56, 793-811. 
FOFONOFF, N.P. 2001 Thermal stability of the world ocean ther- 

moclines. / Phys. Oceanogr. 31, 2169-2177. 
Gnanadesikan, a.. Slater, R.D., Swathi, P.S., & Vallis, 

G.K. 2005 The energetics of the ocean heat transport. J. Climate, 

18, 2604-2616. 

Gregg, M. C. 1987 Diapycnal mixing in the thermocline: a re- 
view. / Geophys. Res. 92-C5, 5249-5286. 

Ledwell, J. R., Watson, A. J. & Law, C.S. 1998 Mixing of 
a tracer in the pycnocline. J. Geophys. Res. Oceans 103, 21499- 
21519. 

Lilly, D. K., Waco, D. E., & Adelfang, S. 1. 1974 Strato- 
spheric mixing estimated from high-altitude turbulence measure- 
ments. /. Appl. Meteor., 13, 488-493. 

LiNDBORG, E. & Brethouwer, G. 2008 Vertical dispersion by 
stratified turbulence. J. Fluid Mech. 614, 303-314. 

Linden, P. F. 1979 Mixing in stratified fluids. Geophys. Astrophys. 
Fluid Dyn. 13, 3-23. 

LORENZ, E. N. 1955 Available potential energy and the mainte- 
nance of the general circulation. Tellus 7, 157-167. 

McEwan, a. D. 1983a The kinematics of stratified mixing 
through internal wavebreaking. / Fluid Mech. 128 41-51 . 

McEwan, a. D. 1983b Internal mixing in stratified fluids. J. Fluid 
Mech. 128 59-80. 

MUNK, W. 1966 Abyssal recipes Deep-Sea Res. 13, 207-230. 

MUNK, W. & WUNSCH, C. 1998 Abyssal recipes II: energetics of 
tidal and wind mixing. Deep-Sea Res. 45, 1977-2010. 

Oakey, N. S. 1982 Determination of the rate of dissipation of tur- 
bulent energy from simultaneous temperature and velocity shear 
microstructure measurements. J. Phys. Oceanogr. 12, 256-211 . 

OSBORN, T. R. AND Cox, C. S. 1972 Oceanic fine structure. Geo- 
phys. Astr. Fluid Dyn. 3, 321-345. 

OSBORN, T. R. 1980 Estimates of the local rate of vertical diffusion 
from dissipation measurements. J. Phys. Oceanogr 10, 83-89. 

Paparella, F. & Young, W. R. 2004 Horizontal convection is 
non turbulent. J. Fluid Mech. 466, 205-214. 

Peltier, W. R. & Caulfield, C. P. 2003 Mixing efficiency in 
stratified shear flows. Annu. Rev. Fluid Mech. 35, 135-167. 

Staquet, C. 2000 Mixing in a stably stratified shear layer: two- 
and three-dimensional numerical experiments. Fluid Dyn. Res. 
27, 367-404. 

Tailleux, R. 2009 On the energetics of stratified turbulent mix- 
ing, irreversible thermodynamics, Boussinesq models, and the 
ocean heat engine controversy. / Fluid Mech., in press. 

Weinstock, J. 1978 Vertical turbulent diffusion in a stably strati- 
fied fluid. J. Atmos. ScL, 62, 3177-3180. 

Winters, K. B., Lombard, P. N., and Riley, J. J., & 
D'Asaro, E. a. (1995) Available potential energy and mixing 
in density-stratified fluids. J. Fluid Mech. 289, 115-228. 

WUNSCH, C. & Ferrari, R.2004 Vertical mixing, energy, and the 
general circulation of the oceans. Ann. Rev. of Fluid Mech., 36, 
281-314. 

f 
f 
f 



12 R. Tailleux: Is mixing efficiency affected by nonlinear equation of state? 



A) New view of energetics of turbulent mixing (Tailleux, 2009) 
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B) Classical view of energetics of turbulent mixing (Winters et al, 1995) 
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Fig. 1. A) New view of tiie energetics of freely decaying turbulent stratified mixing as proposed by Tailleux (2009) versus B) the earlier 
interpretation proposed by Winters et al. (1995). In the new view, internal energy IE is subdivided into a dead part IEq and exergy part 
lEexergy The doublc arrow linking lEexergy and GPEr means that both Wr, laminar and Wr,turbuient can be either positive or negative 
in general. 
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Fig. 2. (Top panel) The squared buoyancy frequency N common to all stratifications considered. (Bottom panel) The thermodynamic 
efficiency-like quantity aP/{pCp) corresponding to the 27 different cases considered. Note that the Fofonoff regime, i.e., the case for 
which GPE decreases as the result of mixing, is expected whenever the latter quantity decreases for increasing pressure. The classical case 
considered by the literature, i.e., the case for which GPE increases as the result of mixing corresponds to the case where the latter quantity 
increases with increasing pressure on average (see Table 1 for more details). 
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Fig. 3. Distribution of the top-bottom temperature difference T{Pm.in) — T{Pmax) as a function of the experiment number. 
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Fig. 4. (Left panel) The dissipation rate of APE as a function of APE, each point corresponding to one particular experiment. Note that 
there is no obvious dependence on the stratification. (Bottom panel) The rate of change Wr,turbuient as a function of APE. This time, each 
stratification is associated with a different curve. 
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Fig. 5. (a) The averaged ratio ^ = Wr,turbuient/ D{APE) as a function of tire experiment number; (b) Tire averaged ratio AGPE / APE 
as a function of the experiment number; (c) The minimum value of d/dz[aP/{pCp) as a function of the experiment number; (d) The 
top-bottom difference of aP/ (pCp) as a function of the experiment number. 
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Table 1. Averaged values of the two ratios ^ = Wr,turbuient/D{APE) and AGPE/APE for the 27 different types of stratifications 
considered in this paper. The quantities [dT /dz\min and AT refer to the minimum value of the vertical derivative of T = aP/{pCp) and 
top-bottom difference of T respectively. S is the salinity used in the equation of state for seawater, T is the mean temperature of the profile 
considered, and Pmin denotes the minimum value of the vertical pressure profile. The top-bottom temperature differences are displayed in 
Fig. 3, while the pressure interval is 10 dbar in all cases. The tabulated values demonstrate that increasingly negative values of ^ coincide 
with increasingly large positive values of dT /dz, as well as with with the increasing importance of non-Boussinesq compressible effects 
associated with an increasing discrepancy between AGPE and APE. The standard case for which f « 1 is achieved close to atmospheric 
pressure. The maximimi negative value of ^ occurs for the lowest S, lowest T, and largest Pmin values considered. 
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